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Henry  Stoomel 
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Introduction 

The  traditional  realm  of  physical  oceanography  has  been  the 
open  ocean,  because,  although  the  observations  even  to-day  are 
sparsely  scattered  in  position  and  season,  the  physical  properties 
of  the  deep  water  seem  to  be  so  steady  that  it  is  possible  to 
make  meaningful  interpretations  by  a  combination  of  observations 
from  many  different,  separate  oceanographic  stations.  Studies  of 
the  oceanography  of  shallow  coastal  waters  is  complicated  by  im¬ 
portant  seasonal  fluctuations  of  properties,  by  tidal  currents, 
by  local  wind  currents,  by  variations  of  fresh  water  run-off  from 
the  coast,  and  by  bottom  friction.  However,  these  shallow  water 
regions  are  of  much  more  economic  importance  than  the  deep  sea 
(Iselln,  1940),  and  many  regional  studies  have  been  made  as  a  part 
of  fisheries  or  pollution  investigations. 

System  of  classification 

The  variety  of  hydrographic,  topographic  and  climatic  con¬ 
ditions  of  coast  lines  is  so  bewildering  that  sweeping  generali¬ 
zations  in  shallow  water  oceanography  have  not  been  made;  shallow 
water  studies  have  been  treated  as  purely  local  phenomena.  Never¬ 
theless,  when  one  reviews  these  numerous  studies  he  is  struck  by 
the  fact  that  there  are  certain  similarities  among  certain  types 
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of  coastal  regions  and  It  Is  natural  to  attempt  to  formulate 
these  similarities  Into  schemes  of  classification. 

The  most  obvious  category  la  that  of  geometric  form: 
broadly  speaking,  one  may  distinguish  between  open  coasts  on 
the  one  hand,  and  small  gulfs,  bays,  lnlets--which  we  may 
speak  of  collectively  as  estuaries — on  the  other  hand.  This 
latter  class  is  the  one  to  which  attention  is  directed  in  this 
article,  and  there  is  a  wide  variety  of  geometric  form,  of  size 
and  depth,  amongst  estuaries. 

But  certainly  the  geometric  form  of  an  estuary  itself  is 
not  of  primary  oceanographic  interest  except  insofar  as  the 
geometry  influences  the  oceanographic  conditions.  For  example, 
it  is  perfectly  possible  to  find  an  inlet  of  only  one  mile  in 
length  which  exhibits  a  salinity  distribution  similar  to  that 
of  an  inlet  100  times  as  long.  From  an  oceanographic  point  of 
view,  therefore,  the  shape  and  dimensions  of  8ui  estuary  have  an 
indirect  influence  only,  and  do  not  serve  as  a  primary  criterion 
for  a  rational  classification  scheme  of  estuaries. 

A  more  rational  scheme  of  classification  might  begin  by  re¬ 
ferring  to  the  predominant  physical  causes  of  the  movement  and 
mixing  of  water  in  the  estuary.  These  may  be  either  tidal, 
meteorological  (wind),  or  river  flow.  For  example,  the  primary 
cause  of  movement  and  mixing  in  the  Raritan  River  (New  Jersey, 
U.S.A. }  seems  to  be  the  tides  upon  which  is  superimposed  a  weak 
river  flow;  in  Pamlico  Sound  (North  Carolina,  U.S.A. ),  on  the 
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other  hand,  the  movement  and  mixing  of  the  waters  seem  to  depend 
entirely  on  the  wind;  and  finally,  the  position  of  the  wedge  of 
salt  water  which  Intrudes  from  the  Gulf  of  Mexico  up  Into  the 
Mississippi  River  Is  controlled  entirely  by  the  volume  of  river 
discharge. 

In  each  of  these  three  cases  a  different  physical  cause 
for  the  observed  distribution  of  oceanographic  properties  pre¬ 
dominates,  but  there  are  estuaries  In  which  It  Is  difficult  to 
discover  a  single  predominating  cause  of  movement  and  mixing — 
for  example,  the  shallow  bays  along  the  northern  coast  of  the  Gulf 
of  Mexico,  where  wind  and  tide  both  play  Important  roles. 

The  tidal  type  of  estuary  has  been  most  studied  from  a  theo¬ 
retical  point  of  view  because  most  harbors  seem  to  be  of  this  type. 
Most  tidal  estuaries  are  fresher  than  the  adjacent  ocean  water, 
but  some,  due  to  small  river  flow  and  large  evaporation,  are  ac¬ 
tually  more  saline  than  the  outside  ocean  water  (e.g. ,  San  Diego 
Harbor,  California,  U.S.A. ).  The  main  effort  of  study  has  been 
directed  at  fresh-water  estuaries,  where  the  velocities  due  to 
river  flow  are  small  as  compared  to  the  tidal  velocities.  Per¬ 
haps  the  most  striking  difference  between  various  tidal  estuaries 
is  the  degree  of  vertical  stratification.  At  one  pole  are  shallow, 
vertically  homogeneous  estuaries  (Raritan  River,  10  feet  deep), 
and  at  the  other  are  deep,  f lord-like  stratified  estuaries  In 
which  the  deep  water  Is  practically  oceanic,  and  only  a  thin  sur¬ 
face  stratum  Is  freshened  (Albernl  Inlet,  Vancouver  Island,  Canada, 
600  feet  deep). 
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Theories  cl  <  tidal  mixing  and  salinity  distribution 

The  flrsTst  rough  approximate  method  of  determining  the  sal¬ 
inity  of  a ti tidal  estuary  from  a  knowledge  of  the  tides  and  river 
flow  la  tie*  "tidal  prism  method" ,  long  In  use  by  engineers  (Met¬ 
calf  and  Idifly,  1935)  In  harbor  studies. 

Denote'  the  hlgh-tlde  volume  of  the  estuary  by  Vg,  the  low 
tide  volunel  by  Vl,  the  volume  of  river  flow  per  tidal  cycle  by 
R.  The  tliaial  prism  la  defined  by 

P  =  Vh  -  Vl 

If  thi  river  water  la  fresh,  and  the  salinity  of  the  ocean 
la  O**.  ties  salinity  of  the  harbor  is  obtained  by  the  assumption 
that  on  each  m  tidal  cycle  a  volume  P  of  water  of  salinity  0^ and 
a  volume  Eotof  fresh,  water,  mix  thoroughly  in  the  harbor,  and  that 
a  volume  1+  -s-R  of  the  mixture  is  expelled.  The  salinity  of  the 
estuary  istidtherefore  P  0^/(P  +  R)  In  the  steady  state. 

KetolB  icn  (1951)  has  Improved  the  simple  tidal  prism  method 
for  computlajing  the  salinity  of  an  estuary,  by  dividing  the  estu¬ 
ary  into  aatmimtier  of  segtnents,  because  It  is  unreasonable  to 
suppose  tht  dt  oomplete  mixing  over  the  entire  estuary  can  occur 
on  each  tile.ss,  It  Is  more  reasonable  to  suppose  that  where  the 
river  flow lihs  weak,  the  maximum  distance  along  the  axis  over 
which  mlxiii;  ss  can  be  oomplete  Is  the  total  displacement  of  a  fluid 
particle  Ij'  ttie  tide.  Using  this  distance  as  the  length  of  his 
segments,  Ss'ntclnum  obtained  remarkably  good  agreement  with  observed 
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salinity  distributions  in  some  shallow,  unstratified  tidal  es¬ 
tuaries. 

Arons  and  Stommel  (1951)  have  translated  Ketchum's  work 
into  the  language  of  physics  of  continue.  The  ordinary  diffusion 
equation  was  employed,  including  a  convective  term  for  the  river 
flow  and  an  eddy  diffusivity.  A,  with  the  tidal  current  velocity, 
Uq,  and  tidal  displacement  ^  q  as  characteristic  velocity  and 
length  respectively. 

A  =  2BUo  ^  o 

where  B  is  a  constant  of  proportionality. 

For  an  estuary  of  rectangular  cross-section,  the  family  of 
curves  shown  in  Figure  1  was  obtained.  Observed  points  for  the 
two  tidal  estuaries  are  shown  on  the  Figure.  The  abscissa  A 
is  the  ratio  of  the  distance  from  the  head  of  the  estuary  to 
its  total  length;  the  ordinate  is  the  ratio  of  the  salinity,  S, 
to  the  ocean  salinity  0  .  The  quantity  F  has  been  called  the 
Flushing  Number 

F  =  Dh2/2B  ^ 

where  D  is  the  river  discharge,  H  the  estuary  depth,  S  o 
amplitude  of  the  tide, CO  the  angular  frequency  of  the  tide,  and 
V  the  volume  of  the  estuary. 

Because  the  curves  presented  here  were  developed  for  a  very 
much  idealized  situation,  it  is  surprising  and  encouraging  to 
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find  tliat  empirical  data  from  actual  surveys  can  be  plotted  on 
tJae  family  with  such  good  agreement. 

An  attempt  to  calculate  the  proportionality  factor  B  from 
the  data  was  unsuccessful,  the  values  being  of  an  order  of  mag¬ 
nitude  different  for  the  two  cases.  Therefore,  it  appears  that 
although  the  shape  of  the  theoretical  curves  is  in  good  agree¬ 
ment  with  the  observations  at  hand,  an  a  priori  calculation  of 
the  flushing  number  is  not  yet  feasible.  It  also  indicates 
that  Ketchum's  analysis  is  a  limiting  case  of  maximum  possible 
lateral  mixing,  which  need  not  always  be  realized.  Neverthe¬ 
less,  the  flushing  number  may  be  a  convenient  concept  to  char¬ 
acterize  estuaries,  Just  as  the  family  of  curves  themselves  is 
a  convenient,  semi -empirical  expression  of  the  mean  salinity 
dlstributl on. 

Dynamics  of  tidal  estuaries 

The  dynamics  of  estuarine  flow  is  not  very  well  understood, 
even  if  we  exclude  cases  of  large  estuaries  such  as  Chesapeake 
Bay  where  Coriolis'  terms  seem  to  complicate  matters,  (as  shown 
both  by  the  mean  salinity  distribution  and  the  Kelvin-type  tide). 
The  flow  in  an  estuary  such  as  Albernl  Inlet  (Tully,  1949),  cer¬ 
tainly  one  of  the  simplest  and  most  clear-cut  examples,  has  not 
been  quantitatively  explained. 

One  might  begin  by  treating  it  as  a  problem  in  free  convec¬ 
tion,  with  an  eddy  dlffusivity  determined  by  the  tides,  but  in¬ 
dependent  of  the  mean  flow.  A  fixed  discharge  of  buoyant  fresh 
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water  is  introduced  by  the  river  at  the  head  of  the  inlet,  and 
carried  away,  by  a  seaward  directed  pressure  gradient,  to  the 
sea.  All  the  time  it  mixes  with  the  underlying  deep  sea  water 
in  the  inlet,  so  that  the  mass  transport  in  the  fresh  layer  in¬ 
creases  greatly  to  the  seaward.  A  compensating  inflow  of  sea 
water  must  occur  in  the  deep  layer  to  replace  the  sea  water 
lost  through  mixing.  A  theoretical  solution  and  understanding 
of  this  problem  (i.e.  ,  a  computation  of  the  salinity  and  vel¬ 
ocity  fields)  has  not  been  achieved  due  to  the  essential  non¬ 
linearity  (Stommel,  1950)  entering  through  the  convective  forms 
of  the  equation  of  salt  transfer.  A  remarkable  feature  of  the 
flow  in  Alberni  Inlet  is  that  the  top  layer  does  not  deepen  to 
the  seaward,  as  might  be  expected,  due  to  the  Increased  mass 
transport  and  mixing.  Instead,  the  surface  water  layer  accel¬ 
erates,  and  even  appears  to  become  slightly  shallower.  It  seems 
possible,  therefore,  that  the  non-linearities  of  the  problem  may 
also  enter  through  the  inertial  terms  of  the  equation  of  motion, 
as  suggested  by  Rossby  (1951),  thus  accounting  for  accelerating 
of  the  surface  flow. 

A  preliminary  model 

To  pursue  this  idea  quantitatively  the  following  rough 


theoretical  steady  state  model  has  been  constructed.  A  deep 

fiord  is  filled  with  ocean  water  of  densitj^  2 

bottom  to  a  depth  5^  =  ^ 2  (^)»  Floating  on  the  top  of  this 


volume  is  a  thin  sheet  of  braoklsh  water  of  density^  ^  ) , 

with  a  free  surface  at  1  coordinate  system 

is  fixed  so  that  the  origin  lies  in  the  deep  water  at  the  head 
of  the  fiord;  the  ^  -  axis  extends  horizontally  toward  the 
ocean;  and  the  2"  “  axis  is  directed  vertically  upward.  To  ob¬ 
tain  a  simple  set  of  equations  the  following  idealizations  are 
Introduced,  in  fact,  they  are  suggested  by  the  nature  of  what  is 
observed  in  real  fiords: 

(i)  The  interface  ^  2  permits  upward  movement  of 

deep  water  into  the  top  layer  where  it  is  mixed;  but  no  mixing 
downward  is  permitted.  This,^  2  independent  of  ^  ,  but 
as  The  mixing  may  be  due  to  winds,  tidal  currents, 

or  the  shear  developed  at  the  interface,  but  in  this  model  it  is 
regarded  as  fixed  Independently  of  the  mean  flow. 

(li)  The  depth  of  the  top  layer,  »  is  very  much 

less  than  that  of  the  bottom  layer. 

(ill)  Within  the  top  layer  mixing  is  so  strong  vertically 
that  the  density  /iC-X)  ,  the  velocity  are  Independent 

of 

For  convenience  we  denote  the  depth  of  the  top  layer  by 
2)  •  and  the  velocity  of  mixing  of  deep  water  into  the 

upper  layer  as  7/m  •  *^^0  mass  flux  per  unit  width  of  the  top 

layer,  ,  increases  by  the  addition  of  deep  water  at  the 

rate  so  that  by  mass  continuity  the  following  relation 
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must  obtain: 


—  (1) 

If  the  salinity  of  the  two  layers  are  denoted  by  S|  ,  and 
,  respectively,  the  salt  transfer  per  unit  width  in  the  top 
layer,  >  increases  by  addition  of  salt  water  from  below, 

and  by  conservation  of  salt  the  following  relation  is  obtained: 


If  the  density  is  simply  a  linear  function  of  salinity,^ 
equation  (1)  and  (2)  combine  to  yield  the  following  form  of 
relation: 


-2^ 


(3) 


Because  of  its  great  depth,  the  velocities  in  the  deep 
layer  are  small,  and  hence  the  horizontal  pressure  gradient  at 
any  depth  within  the  deep  layer  must  vanish.  The  pressure  ^2 
of  the  deep  water  at  depth  is 

p,  .  ^  [j*,  a.  -  -  »)J 

Since  the  vanishing  horizontal  pressure  gradient  is  Independent  of 
within  the  deep  layer,  it  is  convenient  to  consider  it  at  =  0: 
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hi. 


(4) 


or  more  simply: 

^  15) 

The  flux  of  momentum  in  the  upper  layer,  >  is 

not  added  to  by  the  water  mixed  up  from  below  because  by  hypothe¬ 
sis  this  water  has  a  vanishingly  small  horizontal  velocity.  The 
momentum  flux  does  vary  with  -PC  ,  however,  on  account  of  the  pres¬ 
sure  variations.  If  wo  consider  a  unit  length  of  the  top  layer, 
and  compute  the  pressure  at  any  point  within  it  by  the  hydrostatic 
principle,  p, «  (where  is  atmospheric  pres- 


sure),  the  total  pressure  force  over  the  layer  is 
and  the  total  retarding  force  per  unit  length  is 
or  — aji"’  1  u  111  To  this  retarding  force  the  horizontal 
components  of  pressure  acting  on  the  top  and  bottom  of  the  layer 
per  unit  length  and  width  must  be  added:  P.  at  the  top; 

at  the  bottom.  By  adding  these  various  terms  the 


following  equation  is  obtained: 


By  elimination  of  ^  between  equations  (5)  aad  (6)  the 


(6) 
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following  equation  Is  obtained: 


)]- 


Where  V  = 

Equations  (1),  (3)  and  (7)  are  Independent  relations  be¬ 
tween  the  variables  and 


The  meaning  of  equation  (7)  can  be  explained  easily  If  we 
Introduce 


T  ■ 
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a  function  denoting  the  transport  of  the  upper  layer.  Clearly 


•y 


-  v.-n/T 


where  the  subscript  "q"  Indicates  values  of  ^  and  /”  at 
ikO-  The  equation  (7)  then  becomes 


or  If 


OfC.  [is^  ^  T  J 

/(•«  t;  /2 


Changing  variables 

c/LJ^  s 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 
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or  if  4-'  (13) 


The  nuaerical  values  tabulated  In  Table  I  Indicate  that 
there  Is  a  range  of  ,  for  which  the  depth 

of  the  top  layer  decreases  with  Increase  of  transport.  In  the 
case  of  Alberni  Inlet  this  is  apparently  what  happens.  More¬ 
over,  the  value  of  b  increases  rapidly  after  a  certain  value  of 
the  transport,  so  that  evidently  the  two  layer  model  breaks  down 
where 

This  point  of  break-down  of  the  model  appears  to  be  a  kind 
of  "internal  hydraulic  Jump",  which  we  may  see  as  follows: 

The  velocity,  c,  of  an  internal  wave  at  the  interface  is 
given  by 

From  equation  (13)  one  sees  that  the  velocity  of  the  upper  layer, 
u,  is  equal  to  c  at  the  value  b  =  2.  At  points  further  upstream 

^  CL  * 

It  seems  plausible  that  this  region  of  two  layer  break-down 
coincides  with  the  fresh  water  front  so  often  observed  at  the 
mouth  of  deep  estuaries.  Tully  has  communicated  to  the  writer 
the  fact  that  the  Internal  wave  motions  of  the  open  sea  do  not 
seem  to  penetrate  into  the  two  layer  portion  of  an  estuary. 
Apparently  the  "internal  Jump"  acts  as  a  block  to  such  deep  ocean 
waves  and  prevents  them  from  progressing  into  the  estuary. 
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TABLE  1 


b 

0 

.2 


•  4 
.5 

.6 

.8 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

2.2 

3.0 

10.0 

100.0 


M 

0.5 

0.3 

0.125 

0.0 

-0.25 

-0.50 

-1.00 

-1.75 

-3.00 

-5.50 


-13.  00 


+13.0 

+5.0 

2.4 

2.01 
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The  Mixing  Process 

The  chief  unknown  process  involved  in  this  analysis,  and 
for  that  matter  in  any  study  of  deep  estuaries,  is  the  nature 
of  the  mixing  process  of  deep  water  into  the  surface  layer. 

Tully  has  suggested  that  the  mixing  is  due  to  the  action  of  the 
tides,  but  since  in  deep  water  there  is  but  little  vertical 
shear  in  a  tidal  current,  even  in  stratified  water,  it  seems 
to  the  writer  that  the  vertical  shear  of  the  mean  flow  is  not 
entirely  ruled  out  as  being  responsible  for  what  mixing  occurs. 

It  is  not  clear  that  the  mixing  is  really  akin  to  shallow  turbu¬ 
lence  phenomena  because  the  mixing  appears  to  be  a  one  way  affair: 
mixing  up  of  deep  water  into  the  shallow  fresh  top  layer,  but 
little  mixing  down  of  top  fresh  water  into  the  deeper  layers.  This 
appears  to  be  something  quite  different  from  the  ordinary  notion 
of  random  vertical  movements  producing  a  redistribution  of  mass 
in  the  vertical.  Keulegan  (1949)  for  instance,  has  described 
some  tests  on  the  interfacial  mixing  in  stratified  flows  made  in 
flumes.  When  a  current  of  light  water  is  flowing  over  a  resting 
pool  of  heavy  water  with  a  small  velocity  the  interface  appears 
quite  smooth,  but  when  the  velocity  of  flow  of  the  upper  layer  is 
Increased  beyond  a  certain  value,  waves  appear  on  the  Interface 
between  the  two  liquids  of  different  density.  When  the  flow  is 
further  increased,  crests  of  these  waves  break  off  into  the  upper 
fluid  and  are  rapidly  diffused  within  it,  but  there  is  no  corres¬ 
ponding  transport  of  light  surface  water  into  the  lower  water,  and 
he  has  defined  a  certain  critical  velocity  which  we  can  call  u  ,  at 
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which  the  mixing  first  appears. 

Keulegan  has  also  constructed  a  graph  valid  for  a  wide  range 
of  relative  densities  of  the  two  liquids  to  show  the  amount  of 
mixing  as  a  function  of  the  velocity  of  the  current  of  light  liquid 
and  has  obtained  a  relation  of  the  form  »  C  wr 

where  the  mean  flux  of  deep  water  into  the  upper  layer  across 

a  unit  horizontal  area,  and  C, »  3.r 

In  most  natural  estuaries  the  stratification  is  by  no  means 
as  sharp  as  that  which  occurs  in  Keulegan's  flume  experiments,  so 
that  it  is  something  of  a  wild  extrapolation  to  apply  this  formula 
to  vertical  mixing  in  an  estuary.  The  fact  is,  however,  that  it 
gives  reasonably  good  results  consistent  with  the  actual  observed 
mixing  in  Albernl  Inlet.  If  we  make  use  of  this  relation  in  our 
two-layer  model  wo  see  that  by  reasons  of  continuity 
If  we  restrict  ourselves  to  the  region  in  which  V 
then  C  U,  ■  •i  ,  so  that  the  density  of  the  top  layer 

is  given  by  the  relation  /  •  ■y 

Thus  equation  (2)  reduces  to  the  simple  form 

In  the  range  where  D  remains  virtually  constant,  tl»  horizontal 
distribution  of  density  la  given  by  T  •  £  .  Figure  3 

shows  a  comparison  of  the  horizontal  distribution  of  density  ob¬ 
tained  by  using  Keulegan’s  constant  C  and  a  depth  of  the  mixing 
layer  of  7//0  ^ ,  and  the  actual  horizontal  density  distribu¬ 
tion  as  obtained  by  averaging  the  density  in  the  top  layers  in 
Figure  2. 
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